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ABSTRACT 

We present two epochs of observations of TW Hya from the high-dispersion 
near-IR spectrograph ARIES at the MMT. We detect strong emission from the 
Br7 transition of hydrogen, indicating an accretion rate substantially larger than 
previously estimated using hydrogen line emission. The Br7 line-strength varies 
across our two observed epochs. We also measure circumstellar-to-stellar flux 
ratios (i.e., veilings) that appear close to zero in both epochs. These findings 
suggest that TW Hya experiences episodes of enhanced accretion while the inner 
disk remains largely devoid of dust. We discuss several physical mechanisms that 
may explain these observations. 

Subject headings: stars: planetary systems — stars: pre-main-sequence — stars: 
individual (TW Hya) — techniques:spectroscopic 



Introduction 



TW Hya is a nearby (~ 56 pc), young (~ 10 Myr) star ( jWebb et al.lll999l ) surrounded by 
an accretion disk t hat evinces a large inner hole as judged from the obser yed spectral energy 
distribution (SEP; ICalvet et al.ll2002l ) and mid-IR and sub-mm imaging ( iRatzka et al.ll2007l ; 
Hughes et al.l 120071 ). These data can be modeled with an optically thick disk whose inner 
edge is located ~ 4 AU from the centr al star. While observations of 10 /im silicate emission 



(SitkoetaL 



( Eisner et al 



2000l: lUchida et al.l 120041 ) together with spatially- resolved emission at 2 fim 



2006) suggest the presence of dust grains with sizes less than a few microns 



at stellocentric distances R < 4 AU, this inner disk ma terial appears opti cally thin, and 
has been estimated to constitute less than a lunar mass (ICalvet et al.ll2002l ). Such a small 
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amount of material is compatible with the v ery low near-IR veiling measured previously 
fljohns-Krull fc ValentilboOll : lYang et al.ll2005h . 



Warm gas has also been detected within the optically thin inner region of TW Hya's cir- 
cums tellar disk faerczeg et al.l booi iRettig et al.lbood ; ISalvk et al.lbo07l bood ; IPontoppidan et al 
20081 ; iNajita et al.ll2007l ). In fact, the gas appears to be distributed in a di sk with a similar in- 
clina t ion to the optic ally thick dust disk observed at larger radii (i ~ 3-7°; IPontoppidan et al. 



20081 ; iQi et al.ll2004l ). The inner disk of TW Hya, while optically thin, is clearly not com- 



pletely devoid of matter. 

A variety of accretion signatures have been detected in TW Hya, indicating some flow of 
gaseous material through the inner disk. Observations of Ha emission, interp reted in the con- 
text o f a magnetospheric accretion model, indicated M = 5 x 10~ 10 M Q yr -1 (jMuzerolle et al. 



20001 ) . A shock model fitt ed to UV and optical photometry indicated a consistent value of 
M = 4 x 10~ 10 M Q yr _1 flMuzerolle et al.l boooh . A similar model fitted to a sep arate UV 
spectrophotometric dataset yielded M = 2 x 10 -9 M yr -1 (jHerczeg et al.l 120041 ) . Analy- 
sis of two additional epochs of UV data separated by appro ximately two years indicate a n 
accretion rate of ~ 1.5 x 10 -9 M yr" 1 on both occasions (jHerczeg fc Hillenbrand! 120081 ). 
Still higher accretio n rates were suggest ed by modeling of X-ray emission from TW Hya: 
M ~ 10~ 8 M yr -1 (IKastner et al.ll2002l ). In contrast, a subsequent st udy of X-ray emission 
derived an accretion rate of ~ 10~ n M yr -1 (IStelzer &: Schmittll2004r). M odeling of the Civ 



A1549 A line flux yields M = 4 ± 2 x 10" 8 M yr" 1 flValenti et alJl2000h 



While the large dispersion in accretion rate measurements is due in large part to different 
systematic effects for different techniques, in some cases the same technique used on different 
epochs (generally by different investigators) has yielded different accretion rates. Moreover, 
a study of Na D line variabil ity suggested accretion ra tes varying between ~ 10~ 9 and 10~ 8 
M yr -1 over a 1-yr period jAlencar fc Batalhal liooil ). Thus, it seems that the amount of 
material in the inner disk may fluctuate significantly over time. 

In this paper, we present observations of Br7 emission and near-IR veiling. The two 
epochs of data show significant variations in Br7 emission strength, but not in veiling. We 
use these findings to examine the viability to different mechanisms proposed to maintain the 
optically thin inner region of TW Hya's disk. 



2. Observations and Data Reduction 



We observed TW Hya on UT 17 January 2009 and UT 8 May 2009 with the ARIES 
instrument at the MMT. ARIES receives a corrected beam from the adaptive secondary of 
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the MMT, and performs imaging and spectroscopy from 1 to 5 /im wavelengths (although 
the 3-5 /im detector is not yet implemented). We used the echelle mode of ARIES, with a 
resolving power of R = 30, 000. The slit is 1" x O'.'l, well-matched to the diffraction-limited, 
AO-corrected resolution of the MMT. The small slit allows many orders to be placed on the 
detector simultaneously, and our setup included > 25 non-contiguous orders across the H 
and K bands. 

Based on the echelle dispersion, we expect a velocity resolution of approximately 10 
km s _1 . Measured FWHMs of argon lamp lines show that the actual instrumental velocity 
resolution of these observations is approximately 14 km s" 1 . 

TW Hya was observed for 10 minutes on 17 January 2009 and for 40 minutes on 8 May 
2009 (both of these are on-source integration times). In each epoch, observations of TW Hya 
were interleaved with observations, at similar airmass, of an A1V telluric calibrator star, HIP 
54682. The data were flat-fielded using spectra of an incandescent lamp, and extracte d, wave- 



length calibrated, and telluric corrected, using IRAF packages (IMassey et a. 



1992). Wave- 



length calibration used telluric absorption lines and the HITRAN linelist ( iRothman et al. 
20051 ). and is accurate to 1-2 km s _1 . 



We also obtained imaging data of TW Hya and the nearb y star TWA 7 during both 
epochs. TWA 7 is a w eak-line T Tauri star fjWebb et al.lll999l ) of similar age to TW Hya 
(INeuhauser et al.ll2000l ). It exhibits no detectable emission above the stellar photosphere in 
the near-IR, although it does have long-wavelength excess emission indicative of a debris disk 
(IMatthews et al.ll2007l ). Due to the low declinations of both TW Hya and TWA 7, we did not 
have time to obtain spectra of TWA 7. However, since both sources were observed at similar 
airmasses and have similar colors, TWA 7 can be used for relative photometric calibration. 
In January, 2009, both objects were observed using the K s filter with integration times of 0.9 
s. In May, 2009, the atmospheric transparency was better, and these objects saturated the 
detector when the K s filter was used, even for the shortest integration times. We therefore 
observed TW Hya and TWA 7 with a narrowband H2 filter, using 5 s integrations, during 
the May observations. 

During our M ay observations, we also obtained spectra of GJ 568 A, an M3 V star 
(IGray et al.ll2003l ). The total integration time was 25 minutes. The spectra were calibrated 
using the same procedure described above for TW Hya. We use GJ 568 below for a rough 
comparison of the veiling observed in TW Hya with that of an unveiled main-sequence star. 
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3. Results and Analysis 
3.1. Photometry 

For each observed epoch, we computed the counts for both TW Hya and TWA 7 within 
an aperture of 40-pixel (1''6) radius. The sky value was determined as the median of pixels 
with radii between 60 and 80; this value was subtracted from each pixel in the photometric 
aperture. Assuming the 2MASS magnitude of TWA 7, mx = 6.90, we derive the K magni- 
tude of TW Hya at both epochs (Tabled]). The two measurements are consistent with each 
other, and with the 2MASS magnitude, = 7.30. 



3.2. Br7 Spectra 



Our Br7 spectra are shown in Figure [TJ We measure equivalent widths (EWs) from 
these spectra by integrating the line flux between 2.1644 and 2.1668 /mi; this range covers 
all of the Br7 emission seen in both epochs. The main uncertainty in our EW determination 
is the continuum level. We estimate this uncertainty by fitting the continuum in different 
spectral regions. We find 12% error bars for our EW measurements. 

The lines measured in both epochs are blue-shifted. In the January data, the blue-shift 
is more pronounced, and a "blue shoulder" is more prominent in the spectrum. Blue-shifts 
and blue-shoulder-profiles are expected for models of infalling gas, since absorption of Br7 
emission by the infalling gas along the line of sight preferentiall y removes red emission (e.g., 
Hartmann et al.lll994j ; iNajita et al1ll996l : iMuzerolle et al.lll998h . 



The emission in both epochs is substantially broader and more asymmetric than ex 
pected from stellar flares, which typically show strong (Ha line-t o-continuum ratios > 10) 
narrow (FWHM < 50 km s -1 ), symmetric emission line cores (e.g-. lHawley fc Pettersenlll991 



Hilton et al.ll2010l ). Similarly, quiescent chromospheric activity produces emis sion that is nar 
rower and more symmetric than emission produced by accretion processes (IWhite fc Basri 
20031 ). Since the difference spectrum for the two epochs is broad and asymmetric-lacking 
a strong, narrow core-the change between spectra is probably not due to chromospheric 
activity. 

Using the broad-band magnitude of TW Hya derived in §3.1 ^ we can convert the observed 
Br7 EWs into line luminosities. These, in t urn, can be converted into accretion luminosities 
using an empirical relationship derived by ( IMuzerolle et al.lll998l ). To convert the accretion 
luminosities into mass accretion rates requires knowledge of the stellar param eters: M ~ 



L&ccR^/GM*. For this calculation, we use the stellar parameters derived by I Webb et al. 
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( 119991 ). assuming a distance of 56 pc: M* = 0.7 M^ , i?* = 1 R Q , and T e g = 4000 K. These 
are the same values used by iMuzerolle et al.l ( 120001 ). enabling a fair comparison with their 
results. Derived line luminosities, accretion luminosities, and accretion rates are listed in 
Table HJ 



3.3. Veiling 



We estimate continuum veiling in the i^-band using Mg and Al lines near 2.11 /im. These 
lines are found in a spectral region relatively free of telluric lines (unlike, for instance, the 
CO ro- vibrational b andheads), and have been used in past stu dies for veiling measurements 
of young stars (e.g., iDoppmann et al.ll2005l ; Eisner et al.l 120071 ). 



The basic method we use to estimate veiling is to compare the observed spectrum to a 
synthetic spectrum for a stellar photospheric model. We first determine oscillator strengths 
for the Mg and Al transitions by comparing synthetic Nextgen spectra for a solar-type star to 
observed solar spectra. Synthetic spectra relevant for TW Hya are then compute d using stel- 



lar p arameters determined from previous high dispersion optical spectroscopy ( lYang et al. 
2005h : T efr = 4126±24 K, logg = 4.84±0.16, [M/H]= -0.11 ±0.13, and vsini = 5.80±0.63 
km s _1 . The previously derived value of vsmi is smaller than the velocity resolution of 
our data. We rotationally broaden the synthetic spectra by vsini = 15 km s^ 1 , which is 
the convolution of the true vsini (~ 5 km s" 1 ) and the instrumental resolution (~ 14 km 
s _1 ). Finally, we add continuum veiling to these rotationally broadened synthetic spectra 
and find the veiling value that provides the best match between observed and synthetic data. 
Fitted veilings are 0.61 ± 0.09 and 0.57 ± 0.04 in January and May, respectively. The two 
measurements are thus consistent within the 1-er statistical errors. 



Doppmann et al.l ( 120051 ) used this procedure to derive veilings from NIRSPEC data. 



They tested the method with observations of unveiled MK standard stars, and found that a 
correction was needed to produce zero veilings for these stars. They suggested that scattered 
light inside the NIRSPEC instrument may have led to an instrumental veiling contribution. 
Since the ARIES spectrograph has not been characterized in this way, we need to test for 
such systematic effe cts. Using our obs erved spectrum of GJ 568, which is an unveiled M3V 
main-sequence star (IGray et al.l 120031 ) . we can constrain any possible instrumental veiling. 
We estimate the veiling for GJ 568 using the procedure described above with T eS = 3500 K, 
logg = 5, and vsini = 15 km s _1 (Figure [3]). We find Tk = 0.47 ±0.04. 



If we assume that the instrumental veiling is independent of spectral type, at least across 
the M3 to K7 range, then we can use the derived veiling for GJ 568 A to correct the veilings 
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inferred for TW Hya. The corrected veilings for TW Hya are 0.14 ± 0.10 and 0.10 ± 0.06 in 
January and May, respectively. The inferred veilings in the two epochs are c onsist ent with 
each other and with with a previous measurement by I Johns-Krull &: Valentil ( 1200 ll ). 



4. Discussion 

4.1. Inner Disk Variability 

The variation in Br7 emission line-strength between our two observed epochs, and be- 
tween previous measurements in the literature, suggests that the rate at which gas is being 
accreted by the central star is changing with time. However, the amount of dusty material 
in the inner disk does not appear to change significantly, based on the inferred veilings. Our 
measured photometry is consistent with the 2MASS magnitude, providing further evidence 
that the amount of matter in the inner disk of TW Hya is not changing substantially with 
time. 



4.2. Disk Clearing Mechanisms 



Several mechan isms for creating the optically thin clearing in the TW Hya disk have 



orbits (Calvet et al 



been p roffered (e.g., iNajita et al.ll20101) . including planets, which may clear gaps about their 



2002h : photoevaporation of inner disk material (e.g., 



Alexander et al. 



2006 ; Owen et al. 201oh ; or grain growth, which would deplete the population of small grains 



that would produce the near-IR emission. The accretion rates measured 



to reconcile with EUV-driven photoevaporation scenarios ([Alexander et al. 



lere a re difficult 
2006]) , alt hough 



Owen et al. 



Calvet et al.ll2002l ; iHughes et al.l 120071 ) seems to argue 



2010J). The 



not necessarily with EUV+X-ray-dri ven photoevaporat i on models ( 
sharp edge of the TW Hya disk (e.g. 
against grain growth. 

While a planetary explanation may be favored, it too is problematic since the massive 
outer disk of TW Hya should cause rapi d migration and d estruction of planet-mass objects 
unless the disk is unusually inviscid (e.g., Eisner et al. 2006h . While a 10 Jupiter-mass planet 
has been claimed to orbit TW Hya at 0.04 AU (ISetiawan et al.ll2008l ). this would not dynam- 
ically effect the outer disk edge at ~ 4 AU. Moreover, follow-up obs ervations demonstrat e 
that the claimed planetary signal is more likely due to starspot noise ( jHuelamo et al.l 120081 ). 



Any of the mechanisms discussed above are compatible with time-variable accretion. 
Time-variable accretion in the outer disk-for example, due to changes in MRI turbulence- 
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can lead to a larger flow of material through the inner disk, even if the flow is impeded by 
a planet or photoevaporation. In fact, the inn er disk accretion rate can vary indpendently 



since it too may be driven by MRI turbulence ( jChiang fc Murray-Clay! 120071 ). 



While our data indicate a changing accretion rate with time, the low veilings observed 
at both epochs indicate a small and nearly constant amount of dust in the optically thin 
region. A constant dust content with time is also compatible with the unchanging i^-band 
magnitude inferred from our data. 

In fact, the lack of optically thick dust in the inner disk-at any observed epoch-is 
surprising if one makes a simple estimate of the disk column density. For this calculation, 
we assume T disk (i?) oc R^ 1 ^ 2 and E disk (_R) oc R~ 3 ^ 2 . We take the disk scale height to be a 
constant fraction of the radius, -f^disk = 0.1R. The pressure in the disk is thus Pdisk(-R) oc 
Pdisk^disk oc i?~ 3 . 

For the accretion rates derived here (~ 10 -9 M Q yr _1 ), and assuming a viscous disk 
with a = 0.01, we derive a disk surface density at R — 1 AU of 



and a disk density of 



^disk 1n -l2 -3 /o\ 

Pdisk = ~ 3 x 10 g cm . (2) 



If we assume Kdust = 10 3 (appropriate for sub- micron-sized dust; see, e.g. JMiyake fc Nakagawa 



19931 ). and take a gas-to-dust ratio of 100, then the vertical optical depth of the disk at 1 
AU is 

T ~ (^j) ^disk^dust ~ 15. (3) 

Even for accretion rates lower by an order of magnitude, such as those inferred in previous 
studied (Section [1]), r > 1 for this simple, viscous disk calculation. 

To explain the lack of dust in the inner disk, we propose three possible scenarios, all 
of which may operate in the TW Hya system. The first is that the density could be over- 
estimated. An a-disk may not be the correct description for the inner disk of TW Hya, or a 
higher value of a may be appropriate. A gas-to-dust ratio > 100 could also help to maintain 
a higher accretion rate with a smaller column of dusty material. However, observations 
suggest that the gas-to-dust ratio in TW Hya-at least as measured in the outer disk-may 



be substantially smaller than 100 (IThi et al.l 120101 ) . 



The second possible explanation is dust-filtering, where the pressure gradient at the 
inner edge of the optically-thick disk leads to super-Keplerian velocities, and hence trapping, 
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Rice et ahlhood ). By preventing dust 



of dust particles in a certain particle-size range (e.g., 
from reaching the inner disk, filtering can maintain an optically thin region. Moreover, even 
if the gas flow through the inner disk varies, filtering may trap (some of) the additional dust 
particles in the flow, preventing variations in observed veiling. 

Finally, radiation pressure from the star can help to drive dust particles out of the inner 
disk. Alone, this mechanism is unlikely to clear out the inner disk or TW Hya, since an 
optic ally thick disk can be replenished by accretion faster than material can be blown out 
(e.g., iTakeuchi fc Lin! 120031 ). However, if the inner disk is optically thin to begin with, for 
example due to dust filtering, then radiation pr essure acts directly on all dust i n the system 



and c an be an efficient means of removal (e.g., ITakeuchi &: Artymowiczl l200ll ; Eisner et al. 



20061 ). 



Radiation pressure and dust filtering act in a complementary way. Dust filtering is most 
efficient for particles with sizes of tens of microns, since these are not too well coupled to 
the gas but don't have the inertia of larger bodies (IRice et al.ll2006r). Radiation pr essure is 
most effective for smaller particles, with sizes < 1 /im (e.g.. IWeidenschillingl Il977| ). Thus, 
small particles that pass through the "filter" into the inner disk may be pushed back out by 
radiation pressure. The combination of filtering and radiation pressure may be sufficiently 
efficient to maintain an inner disk virtually devoid of dust even during periods of enhanced 
gaseous accretion. 
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Table I. Derived Properties 



Epoch 


m K 


Br 7 EW (A) 


LBr 7 (L©) 


Lace (L©) 


M (M Q yr- 1 ) 


r K 


17 Jan 2009 
8 May 2009 


7.28 ±0.02 

7.29 ±0.05 


-7.5 ±0.9 
-3.6 ±0.4 


(3.4 ± 0.5) x 10" 5 
(1.6 ±0.3) x 10- 5 


0.06 ±0.01 
0.03 ±0.01 


(3.0 ±0.5) x 10" 9 
(1.2 ± 0.3) x 10- 9 


0.14 ±0.10 
0.10 ±0.06 
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Fig. 1. — Spectra of Br7 emission from TW Hya observed with ARIES at the MMT on 
17 January 2009 (solid histogram) and 8 May 2009 (dotted histogram). Telluric absorption 
features have not been divided out from these spectra, and the velocities of telluric lines 
in both epochs are indicated by vertical lines. The stellar velocities used to create the 
velocity axes in these plots were determined by cross-correlating observed spectra of the CO 
bandheads with model spectra, and are accurate to ~ 5 km s _1 . 
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Fig. 2. — Spectra of TW Hya (solid black histograms) and model spectra that have been 
rotationally broadened and veiled (thick gray curves), for 17 January (left) and 8 May (right). 
We also plot the unveiled synthetic spectra with gray dotted curves. As discussed in §3.3[ 
systematic errors, in particular instrumental scattering, lead to inferred veilings larger than 
true values. 
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Fig. 3. — Spectrum of GJ 568 A (solid black histogram) and a model spectrum that has been 
rotationally broadened and veiled (thick gray curve). We also plot the unveiled synthetic 
spectrum with a gray dotted curve. 



